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a multielectron target as OCS. 
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I. INTRODUCTION 

The carbon oxysulfide (OCS) molecule is well known to be of considerable importance for the 
role it plays in the carbon chemistry and sulphur chemistry chains of reactions in molecular as- 
trophysics environments like those in the diffuse and dark interstellar clouds (DISC)[1]. It is also 
of interest in the realm of cold molecular plasmas because of its function as a reaction quencher 
into the plasma formation. It follows, therefore, that collisional processes induced by low-energy 
electrons in this molecular gas have triggered the interest of both experimentalists and theoreti- 
cians, intrigued by the possible role that the permanent dipole of this target molecule can play in 
characterizing the scattering attributes. 

On the theoretical side, earlier calculations were carried out by Lynch et al.[|2D using their con- 
tinuum scattering model, with which they found the presence of a II symmetry shape resonance 
around 1.15 eV and a sharp increase at threshold of the £ symmetry partial cross section. The dif- 
ferential cross sections (DCSs) for elastic scattering and for vibrational excitations were measured 
by Sohn and coworkers |3|], where the authors also observed a n* resonance and a strong excita- 
tion of the bending mode associated to dipole coupling and to a possible bound state of (OCS~)*. 
Further measurements are due to Szmytkowski et al. i4\ who measured the total, integral cross 
sections of the title molecule and confirmed a il-type resonance around 1.15 eV. More recent mea- 
surements have been carried out by the Japanese group [50, who compared similar measurements 



with electrons and positrons as projectiles and confirmed the presence of the %* resonance around 
1.15 eV. The vibrational excitation of OCS has been measured by another study from the same 
group [60 . Calculations of elastic cross sections, integral and differential, and of the momentum 
transfer cross sections have been performed by Bettega et al. [71], who employed the Schwinger 
multichannel method with pseudopotentials within the static-exchange approximation and exam- 
ined a range of energies from 5 to 50 eV. 

A combined theoretical-experimental study has been conducted by Michelin et al. [8] who car- 
ried out both measurements and calculations in the energy regime from 0.4 eV to 600 eV. Further 
calculations of both integral and differential cross sections were reported very recently by Bettega 
et al. [9], who employed again the Schwinger multichannel method with pseudopotentials,adding 
the polarization interaction to their static-exchange initial scheme. A further "addendum" to the 
above results was recently published with regard to the S-symmetry partial cross sections Jioll . 

In the present study we intend to discuss in some detail the behavior of the computed angular 
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distributions with respect to the available experiments over a very broad range of collision energies, 
using an entirely different method as that employed by the earlier, just as extensive, calculations 



of references 



The following Section describes briefly our theoretical and computational 
approach while Section III reports our results for both the partial integral cross sections and the 
differential cross sections. The last Section IV summarizes our present conclusions. 

II. THE THEORETICAL MODELLING 

A. The separable exchange representation 

One of the major requests from accurate calculations of low energy electron-molecule scatter- 
ing observables is the correct description of the non-local part of the interaction potential. For 
electronically elastic collisions, a physically realistic treatment usually starts with the static in- 
teraction and views the exchange potential as a correction to the latter due to the imposition of 
anti-symmetry on the total electronic wavefunction that describes the scattered electron plus the 
target bound electrons. Even in this case, however, the correct treatment of non-local forces for 
polyatomic molecules still represents a rather formidable task and therefore increasing use is be- 
ing made of separate representations which expand the exchange kernel in terms of products of 



one-electron functions Jilt]. This approach will be adopted in the present paper, following our 
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extension of this method to polyatomic molecular targets [|12Lll3|l . It has the advantage of being a 
non-empirical method for treating exchange forces while still offering computational savings with 
respect to the exact solution of the integro-differential equations of non-local electron-molecular 



scattering theory Ill4ll . The earlier work on separate exchange in electron collision with small 



diatomics has been rather encouraging since fairly small exchange basis sets added little more 



121] and therefore we use such an approach 



computer time than a local potential calculation |jllL 
on the linear polyatomic target of the present study. 

One starts by approximating the exchange potential W(r, r e ) by the truncated separate form 

N 

W(r,r e )^5> a (r)W Qi/3X/3 (r e ) (1) 

a,f3 

where the {x a } now constitute an additional, new set of Cartesian GTOs not orthogonal to each 
other, nor to the occupied molecular orbitals (MOs) of the target SCF basis set which describes the 
ground state of the target molecule. The vector coordinates (r) and (r e ) describe the positions of 
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the bound and of the scattered electrons, respectively. The W aj p now represents the two-electron 
interaction operator over the truncated, discrete basis. 

When performing the separable expansion (1) care should be taken of how to specify the ex- 
change basis functions {xa} in order to avoid taking too large an expansion so that the calculations 
become too massive or begin to suffer from linear dependency effects. We have carefully checked 
this aspect of the problem by analysing in each case the corresponding eigenvectors of the overlap 
matrix and modifying the basis set accordingly. 

The required exchange matrix elements for the bound MOs of the target, taken to be a closed- 
shell structure, are therefore given by first calculating the following matrix elements ,via a standard 
bielectronic subroutine 

B 7T = J dr J dr eV ? 7 (r)^(r,r e )v? r (r e ) (2) 

In the above equation {v? 7 } denotes the set of doubly occupied self-consistent- field (SCF) target 
MOs. Hence, using (1) one can further write 

N f f 
B 7r = dT(p y (T)xa(r)W at p J dr e p T (r e )x /3 (r e ) (3) 

a,/3 

our further compact the overlap integrals 

N 

B7T = ^2 S~faWa,(3Sp T (4) 
a,P 

One can finally obtain the bound-continuum scattering matrix from the separable representation 
by writing it down via the following product of the matrices already defined in the above equations: 

W = S^KS" 1 (5) 

where the S ja are the overlap matrix elements between the additional GTO set and the original 
expansions describing the bound MOs. 

B. The scattering equations 

Within a single-centre expansion (SCE) of the continuum wavefunction and of the interaction 
potential, the use of the present static-separable-exchange approximation gives rise to a set of 
coupled integro-differential equations 
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o(r e ) = Y,\ U ii'( r e)wi (r e ) + Y,® l a(re)W a p [ dr^(r)u l% (r)\ (6) 



d 2 1(1 + 1) \ 
-zpi + — — u ^ 



where 



M*e) = | 5r(f e )y(r e )^'(r e )rff e (7) 

and 

K(r e ) = (r e ) / rff 3 ^(r e )5r(r e ) (8) 

which integrates over real spherical harmonics to yield the radial part of each new GTO function. 
Here the ip's are the orbital used in equation (2). 
Furthermore 



Sr(r e ) = ^OTfc) ± ("inr m (r e ) (9) 

with the parity index p=0 or 1 and where the same equation (8) also holds for the r variable in (6). 

One can now express the solution as a linear combination of homogeneous and inhomogeneous 
terms: 



«IIo(»-e)=«S (re)+E«?( r e)Cfe (">) 



where 



{k 2 - H l )ul(r e ) = ]T U w (r e )u vl0 (r e ) (11) 
v 

(k 2 - H l )u?(r e ) = J2U u >(r e )uZ(r e ) + & a (r e ) (12) 
v 

The coefficients Cf are found to satisfy a set of linear equations 

J2A aP C^B al0 (13) 

where 
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and 



A*/3 = <W - W <*i [ ®p(r e )u?,(r e )dr e (14) 



= E ^ /° ^(^X (r e )rfr e (15) 
i'p 

The final numerical integration of the ensuing Volterra equations was then carried out as already 



described in Jain and Norcross [1511 . 



HI. COMPUTED AND MEASURED CROSS SECTIONS 
A. Computational details 

The interatomic distances were fixed at their experimental values of R oc = 2.196 a and R cs = 
2.941ao- The two components of the dipole polarisability tensor were taken also, in the asymptotic 
interaction, to be the experimental values: a (R eq ) = 35.1ajj and a 2 (R eq ) = 17.5ajj. The above 
choice was dictated by our desire to realistically treat low-energy effects which are known to be 
strongly affected by polarisabilities. 

To describe the target electrons in the neutral ground state we employed a double-zeta D95V 
basis set plus polarization orbitals as in the expansion described as D95V* 01. The obtained total 
energy was therefore -1020.4122528 hartrees, with a dipole of 0.32 a.u. and a quadrupole of 0.87 
a.u., to be computed with the experimental values of 0.27 a.u. for the dipole L4]. No experimental 
value of the molecular quadrupole was found to be available. 

The radial grid chosen for the Volterra integrals went from 0.005 a to 10.0 a , beyond which 
the integration was extended out to 1,000 a using the asymptotic static multipoles plus the po- 
larization potential. The multipolar went up to X max =70 while the bound orbitals were expanded 
around the center of mass up to l max =ll&. 

Table I reports the additional GTO's employed to describe the separable exchange and used in 
the dynamical calculations outlined in the previous Section. 

The single-center partialwave expansion of the discretized, continuum orbitals employed to 
treat exchange involved l ma x=l5 for the symmetry, l max =20 for the II symmetry and l ma x=15 
for the A symmetry. 
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TABLE I: Basis set functions (types, locations and exponents) for the continuum orbitals. 



location: on the 


on the C 


on the S 


on the Center of mass 


V: 


' 4 


Y: 


Y: 

C 4 


s:8.0,4.0,2.0, 


s:8.0,4.0,2.0, 


s:8.0,4.0,2.0, 




1.0,0.5,0.25 


1.0,0.5,0.25,0.12 


1.0,0.5,0.25,0.12 




p z :2.0, 1.0,0.5,0.25 


p z :2.0,1.0,0.5,0.25 


p z :2.0, 1.0,0.5,0.25 




d xz : 1.0,0.5,0.25 


d xz : 1.0,0.5,0.25,0. 12 


d xz : 1.0,0.5,0.25,0. 12 




II: 


II: 


IT: 


IT: 


p x :32.0,16.0,8.0, 




p x :32.0,16.0,8.0, 


p x :32.0,16.0,8.0,4.0, 


4.0,2.0,1.0,0.5 




4.0,2.0,1.0,0.5 


2.0,1.0,0.5 


d xz :4.0,2.0, 1.0,0.5,0.25 




p xz :4.0,2.0,1.0,0.5,0.25 


d xz :4.0,2.0, 1.0,0.5,0.25 


A: 


A: 


A: 


A: 


d xx : 16.0,8.0,4.0,2.0, 




d xx : 16.0,8.0,4.0,2.0, 


d xx : 16.0,8.0,4.0, 


1.0,0.5,0.25 




1.0,0.5,0.25 


4.0,2.0,1.0,0.5 


f x 2 z :4.0,2.0, 1.0,0.5,0.25 




f x 2 z :4.0,2.0,1.0,0.5,0.25 


d xz : 4.0,2.0,1.0,0.5,0.25 



The above set of functions was the result of a series of numerical tests at a set of chosen energies 
where we increased the number and quality of the additional GTO's until the corresponding, final 
K-matrix elements varied by less than 0.1%. 



B. The integral cross section 

The results reported by figure 1 analyze the partial and total integral, elastic cross sections 
(rotationally summed) for their overall energy dependence over the broad range of energies which 
is experimentally available |3, 4, 5], but we only report a comparison with some of the data in 
order not to crowd the figure excessively. We also show in the lower panel the energy behavior of 
the individual partial cross sections for the contributing partial symmetries. 

It is reassuring to see that the present calculations follow closely the experimental findings from 
below 1.0 eV of energy and up to 100 eV: this is a rather good result considering the complexity 
of the molecule and the broad range of energies spanned by the calculations. Furthermore, we 
see that the resonance position is obtained reasonably well from our calculations, albeit with a 
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FIG. 1: Upper panel: Computed and measured total elastic integral cross sections as a function of energy. 

Open circles: expt.s from ref.[3]; open triangles: 
expt.s from ref.[5]. Lower panel: partial cross sections for the S-(solid line), Il-(dashes) and A-(dots) 
symmetries. 



Solid line: present results: Dots: calculations from ref. 



sharper peak due to the absence of vibrational averaging, which also come very close to the recent 
calculations of Bettega et al. [9] that employed an entirely different method for obtaining it. 

The lower panel in the figure reports the energy dependence of the various partial cross sections 
which contribute to the scattering process. The calculations clearly show that: 

1 . the resonance around 1 eV is due to the II symmetry i.e. to the expected dynamical (cen- 
trifugal) trapping of the electron into a ir* metastable orbital by the /=1 angular momentum 
barrier; 

2. there is a second n* resonance at higher energies (~ 10 eV) which is a broader one and which 
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is also seen in the experimental data; 

3. the S-symmetry exhibits the expected dominance of s-wave scattering as one moves down to 
low energies [10] and the corresponding cross section goes through a Ramsauer-Townsend 
minimum around 0.7 eV, as also analyzed and discussed by recent calculations of Bettega et 

al. H]; 

4. the E-symmetry cross section also suggests the presence of a a* resonance around 4-5 eV 
which is also seen in the experiments (upper panel) as a shoulder on the raising cross section 
energy dependence around 4 eV. In Fig. 2 we report the behaviour of the eigenphase sum 
associated with the S symmetry component. We see that, at low energies, this quantity goes 
through zero, as in ref. Jlp| . and therefore it suggests the presence of an RT minimum in 
the elastic cross section. It further increases from zero energy values and goes through a 



maximum as shown by similar data in ref ill On 



In conclusions, the present calculations for partial and total integral cross sections indicate 
that our present approach can reliably describe the experimental behavior over a very broad 
range of energies and can also reproduce earlier calculations very closely to their findings 
as far as partial symmetry contributions are concerned. 

The previous literature flflfl 

also discusses the possibility of having a virtual state close 
to zero energies detected by a negative value of the scattering length. Since the present 
system is a polar target, no modified effective range theory can be applied jlq] and no 
scattering length could therefore be uniquely defined. Different procedures therefore need 



to be used to extract information on the possible existence of virtual states LI 8] in polar 
targets. In the present case, however, we decided that such a study was outside the scope of 
our current aims. 



C. Angular distributions at low energies 



We have shown above that the well marked, low energy resonance appears, but in experi- 
ments and calculations, around 1 eV of collision energy, and therefore it becomes of interest 
to also see what the angular distributions look like across that range of energies. The results 
of our calculations, and their comparisons with the existing experiments, are shown in the 
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various panels of figure 3 and figure 4 where the energies vary from 0.4 eV up to 3.0 eV. 
The experiments are those of Sohn et al. [3] and from Tanaka y5|], respectively marked by 
open circles and open triangles in all the panels, while the calculations (dashed curves) are 
from reference [|80. The following comments could be made: 

(a) at the lowest collision energies the scattering is dominated by the weak dipole interac- 
tion which causes a marked forward scattering behavior, which is also well reproduced 
by our calculations. On the other hand, in the larger angle region the agreement be- 
comes less satisfactory, possibly due to our model treatment of short-range correlation 



effects. We also see that at 1.15 eV the earlier calculations QSD also do not agree well 
with experiments. 

(b) as the collision energy moves across the position of the strong n* resonance, as shown 
by the panels of figures 3 and 4, we see that the angular distributions become much 
flatter from $ ~ 20° and out to 140°, with a markedly reduced forward peak, with the 
latter resuming as soon as the energy leaves the resonance region (e.g. see panels in 
figures 3 and 4). 

(c) We also see from Fig. 4 that, while the agreement with experiments is good at 2.0 eV, 
and remains reasonableat 2.5-3 eV, indicating numerical convergence of our present 
results to the same level of reliability with respect to the experiments. 

In other words, we see the dominant presence of a forward scattering behavior of the DCSs un- 
less we are at the resonance position, where the 'orbiting' features of the trapped electron distorts 
the effects coming from the permanent dipole of the target molecule. 



D. The DCSs at higher energies 

Since the experimental data for the angular distributions are available over a broad range of 
collision energies, we have analyzed them rather carefully from above the dominant n* resonance 
out to about 100 eV. The comparisons with the measured data are shown in figures from 5 through 
7. The following comments could be made from a perusal of the data shown in those figures: 

1. the present calculations are seen to follow measurements remarkably well, both in shape 
and size, over the whole energy range; 
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FIG. 2: Computed eigenphase sum behavior within the S symmetry component of the scattered electron. 



2. at 5.0 eV the data from Sohn et al.[2] differ from those of Tanaka et al. [5] in the small- 
angle region since they indicate there a strong forward peaking of the angular distributions: 
our calculated values follow those measurements very accurately, thereby confirming the 
experimental findings in that angular regime; 

3. as the collision energy increases one sees an increasing flattening of the DCSs as a function 
of the scattering angle and the appearance of the forward peaks at increasingly smaller 
angles: the calculations follow suit in the sense that indeed show the same general behavior 
as that indicated by the experiments. 



IV. PRESENT CONCLUSIONS 

The calculations discussed in this work have analyzed, using a quantum treatment of the scat- 
tering dynamics, the elastic cross sections for electron impact on the OCS target molecule in the 
gas phase. In particular, we have employed the exact description of the static potential and a sep- 
arable representation of the exchange potential, both represented via a single-center expansion, 
and we have further added correlation-polarization effects via a density functional formulation 



that we have discussed many times before II 1711 . The corresponding integro-differential equations 
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have been solved via quadratures of Volterra equations as discussed in section II and the angular 
distributions have been computed including the Born dipole corrections beyond the partialwave 
value of / max =30 Q20]. The final results for the integral cross sections indeed confirm the presence, 
around 1 eV of collision energy, of a narrow and intense shape resonance of Il-symmetry associ- 
ated to a well-known 7r* resonance for the title system (e.g- see discussion in ref.sjv, s|]. It also 
suggests the presence of further two resonances of ix* and a* symmetry and larger widths, together 
with a clear Ramsauer-Townsend minimum in the cross sections around 0.7 eV Jioll . Furthermore, 
the various available angular distributions for the elastic scattering have also been analyzed over 
a very broad range of collision energies, spanning nearly 100 eV, and compared with the existing 
experiments and with the earlier calculations [sHaO]- The comparison of all the distribution data 
were presented in the previous Section and they show fair agreement between the measurements 

nun 

and with other computations [EISLllOD- The present calculations also show particularly good cor- 
respondence between computed and measured distributions for energies from 3.0 eV and up to 
100 eV. Considering the complex, many-electron structure of the target and the broad range of 
collision dynamics which has been analyzed, our results do indicate the robustness of the adopted 
dynamical integrator and the reliability of the theoretical modelling which we have employed in 
this study. 
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FIG. 3: Computed and measured angular distributions. Solid line: present calculations, open triangles: 
expt.s from ref.yfl. Open circles: expt.s from ref. [5], dashes: calculations from ref. [8]. Upper left panel: 



0.4 eV. Upper right panel: 0.6 eV. Lower left panel: 1.15 eV. Lower right panel: 1.5 eV. 
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FIG. 4: Same as in figure 3 but for different collision energies. The dashed lines are calculations from 
ref . Q8I] . Upper left panel: 1.7 eV; upper right panel: 2.0 eV; lower left panel: 2.5 eV; lower right panel:3.0 
eV. 
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FIG. 5: Computed and measured angular distributions. Solid lines: present calculations. Dashes: calcula- 
tions from ref.[8]. Open circles: expt.s from ref.[5]; triangles: expt.s from y|]. Upper left panle: 5.0 eV. 
Upper right panel: 6.0 eV. Lower left panel: 7.0 eV. Lower right panel: 8.0 eV. 
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FIG. 6: Same as in figure 5 but for different collision energies. Upper left panel: 9.0 eV; upper right panel: 
10.0 eV; lower left panel: 15 eV; lower right panel:20.0 eV. 
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FIG. 7: Same as in figure 6 but for different collision energies, Upper left panel: 30.0 eV. Middle panel: 



60.0 eV. Lower panel: 100.0 eV. 



